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Phase separation in systems with absorbing states. 
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Cy ' Abstract. - We study tiie problem of phase separation in systems with a positive definite 

order parameter, and in particular, in systems with absorbing states. Owing to the presence 
of a single minimum in the free energy driving the relaxation kinetics, there are some basic 
• properties differing from standard phase separation. We study analytically and numerically 

(~| ' this class of systems; in particular we determine the phase diagram, the growth laws in one and 

Q ^ two dimensions and the presence of scale invariance. Some applications are also discussed. 
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The presence of conservation laws plays a key role in determining the phenomenology and 
critical behavior of phase transitions both in equilibrium and nonequilibrium systems. For 
instance, the dynamics of models for magnetism change dramatically depending on whether 
the order parameter is a conserved quantity or not ||l[. In the conserving case, models may 
exhibit phenomena such as spinodal decomposition (SD hereafter) and droplet growth via 
condensation-evaporation H, ^, ^, which are absent in models with non-conserving dynamics. 
The problem of SD that appears, for instance in binary alloys, polymers and spin systems 
with conserved magnetization, have attracted a lot of interest in the last decades ^, || , and 
despite the efforts devoted to understand their phenomenology, some of their aspects remain 
unclear. The key issue is that of understanding the dynamics of a system quenched from a 
state in the homogeneous phase into a broken-symmetry phase. Depending on the value of the 
conserved order parameter, M, different types of morphologies can show up. The archetypal 
field theory representing this class of systems is the model B (or Cahn-Hilliard equation), 
defined by the following Langevin equation 



,26Fmx,t)] 



dt<j,{x,t)^^' '2, +V-7?(x,t) (1) 
d(p[x, t) 
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where F[(j3(x,t)] is a free energy functional given by 



F[(j)] = j dx 



(2) 



with a and h positive constants, and ri{x,t) is a d-dimensional noise vector whose components 
obey 

{r],ix,t)) =0, {T^,{x,t)7j,ix',t')) ^ D6^^Hx-x')6{t-t')S,,,. (3) 

A brief discussion of its phenomenology is as follows: an homogeneous initial condition with 
averaged field density (j)Q, located in the unstable zone of the free energy, experiences a linear 
instability which favors the growing of Fourier modes with fc 7^ (waves), and the initially 
homogeneous phase roughens, in such a way that locally the system tries to relax to one of the 
two minima of the free energy, conserving the total order parameter. Let us point out that the 
presence of two different minima in Eq. (||) is a key ingredient for phase separation. After an 
initial stage the wave amplitude grows, until it saturates, i.e. the maxima and minima reach 
the maximum and minimum of the free energy potential, and the symmetry with respect to 
the initial value of the field is lost Afterwards the clusters of the minority phase coalesce 
via diffusive annihilation of interfaces. For metastable initial conditions the system is locally 
stable, but fluctuations can take the local order parameter to one of the minima, in this way 
droplets are formed, and they grow via nucleation. 

In a seminal work Lifshitz and Slyozov and independently Wagner |Q, derived some 
exact results for very diluted quenches, i.e. when the volume fraction of one phase is much 
larger than the other one. For example, the averaged gyration radius of clusters in the minority 
phase, i?, grows with time as i? « t^^'^, while it grows like i^/^ in the non-conserved case. These 
results have been shown to be independent of the spatial dimension d ^ and confirmed in 
numerical studies Although the Lifshitz-Slyozov- Wagner theory (LSW) is only valid 

for very dilute systems, it has been argued that a similar law t^^^ is valid also for generic 
volume fractions . An intermediate regime has also been identified in some cases ; it is 
characterized by the presence of soft domain walls and a growth law t^/^; after this regime the 
domain walls harden (become sharper) and the real asymptotic behavior t^^^ sets in. In other 
cases in which bulk diffusion is stopped the assymptotic domain growth law is characterized 
by an exponent 1/4 Q. 

Another important result is the identification of a single relevant length scale in the late 
stages of phase separation. It is given by the averaged droplet size. This makes it possible to 
describe the late stages of growth in terms of scaling functions [||, p^, and makes the 

problem suitable to be studied using renormalization group techniques |16~Q. It is accepted 



that in the scaling regime stochastic effects are irrelevant and therefore the same behavior is 
obtained in purely deterministic systems with the same symmetries and conservation laws. 

In this paper we face the generalization of the previous ideas to a vast class of nonequilibrium 
systems, namely, the class of systems with absorbing states. It is well established that all the 
systems exhibiting a continuous transition into a unique absorbing state, without any other 
extra symmetry or conservation law, belong to the same universality class, namely, that of di- 
rected percolation ^ . Among other models in that broad class are the following: directed 
percolation, the contact process, catalytic reactions on surfaces, branching annihilating random 
walks, damage spreading and self organized criticality (see |l^ and references therein). The 
Reggeon field theory (RFT) is the minimal continuous theory in this vast universality class 
[ p!8| . In its conserved version it is defined by a Langevin equation analogous to Eq. (^ where 
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now F[(j){x, t)] is a different free energy functional given by 




(4) 



and 



D(f?{x,t)S^'^\: 



X - x')S(t 



(5) 



The main properties of this equation in contrast to those of model B are: (i) all the dynamics 
ceases for (j){x) = 0, as corresponds to a system with an absorbing state; (ii) the RFT describes 
the dynamics of a density like, positive definite, order parameter, and (iii) no up-down 
symmetry is present in this type of systems and therefore the first nonlinearity in the free 
energy is a cubic term instead of a quartic one as in model B (observe, however, that none of 
the forthcoming results depend on the degree of the non-linearity). 

Owing to these last two differences it is clear that the essential mechanisms present in 
standard phase separation cannot straightforwardly be translated to the conserved RFT: The 
deterministic free energy F exhibits only one minimum (at if a < 0, and at = a/b if a > 0), 
and the process of phase separation, if present, does not correspond to a local relaxation of 
the density to the two minima of the local free energy as in model B. 

Before proceeding further, let us present a simple toy model in economics that shares the 
properties of our previously described system (as described by Eq. (|l|), (^) and (||). This 
toy model can be view as a discretization of our proposed equations, and maybe be useful to 
develop an intuition of the forthcoming conclusions and results. The model is defined on a 
2-dimensional discrete lattice, each site representing an individual. A field variable describes 
the wealth of each individual at each time. The wealth is updated at each time step following 
three different mechanisms (1) the wealth is driven by an external agent (a fair government) 
that tries to keep all the individual equally rich and corresponds to the potential in our model; 
(2) the wealth is redistributed periodically among nearest neighbors ( this corresponds to the 
Laplacian in (^ and describes fair citizens), and (3) in addition there are random fluctuations. 
The dynamics is local (wealth flows only through nearest neighbors) and conserves the total 
wealth. The problem is whether the individuals will remain equally rich on average or not 
with these Utopian rules. As we will show in what follows the answer depends on the total 
wealth of the society. 

Coming back to the general problem, a straightforward analysis of local stability of our 
continuous model shows that the homogeneous state is unstable with respect to fluctuations 
whenever </> < a/2b. The problem arises of what is the kinetics in this case; if the system 
segregates in different phases what do they correspond to? 

By Fourier transforming the Langevin equation, it is easy to verify that the fastest growing 
mode is km — y/ {a — 2b(po)/2, and therefore the fastest growing length, given by A = 2n/k„i 
is A = 2-K I ^{a — 260o)/2. The system size has to be larger than A for the instability to be 
observed; for flnite systems of size L, the effective instability point is shifted from (/)o — a /2b 
to a smaller value ~ {a —j-^)/2b, determined by the condition fc„, = L. 

For parameter values in the locally unstable regime initial homogeneous conditions break 
down, and a wavy solution dominated by the wave vector k.„i sets in. This wavy solution cannot 
grow in amplitude indeflnitely due to the physical constraint at = 0; and, in contrast to what 
happens in standard SD the maxima and minima of the wavy solution do not collapse to the 
two minima of the deterministic potential (which in our case is monostable). Therefore, if the 
system is unstable the question arises of what is the true stable solution. In order to answer 
the same question for the model B, one can solve the equation 5F[(j)\/5(j) = with F[(j)\ given 
by Eq. (j^). A well known solution of the previous equation is 0(x) = yja/bian\i{y/ax/ (26)). 
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This is a topological solution connecting the two minima of F. In our case such a topological 
solution cannot exist due to the presence of only one minimum in F. This, however, does not 
mean that the full stationary equation, i.e. W'^{SF[<j)]/d(j)) — does not admit a topological 
solution, even though in this case it is a much less intuitive one. Let us look for such a solution, 
and let us do it first in c? = 1 (the generalization to higher dimensions being straightforward) . 
The general solution of the Laplace equation, V^/(a;) — in d — 1 is f (x) — h + gx where 
h and g are undetermined constants; therefore a solution of the stationary problem is also a 
solution of — 1^ = h+gx (where we have just replaced / by —|^). Imposing periodic boundary 

conditions, we can fix 17 = 0, while h is specified by the initial condition, hV = — / dx^\^=^„ 
where V is the system volume (which transforms into h ~ a(j){) — bipQ with (f>o ~ (f){t — 0) 
for homogeneous initial conditions). Observe that the previous equation, — h, can be 

rewritten as ^ = 0, with G[4>\ — F[(j)\ + hcj), and therefore the effect of the initial condition is 
to originate a constant effective field, /i, coupled linearly to the field variable. We can now look 
for solutions of this modified equation, which may have two minima. In fact, the morphology 
of the new deterministic potential (in terms of G) is as follows: 

(i) For (/)o > a/6, h is negative, and the effective potential has only a minimum (there is 
another, nonphysical extremum of the effective potential located at some negative value of (j)). 

(ii) At a/h > (jjQ > a /2b the effective potential develops a maximum at = o-/b — (t>o, 
and the minimum is at (f>m = </>o- Between = and 0m the potential is a monotonously 
increasing function of 0, therefore at = the potential exhibits a relative minimum even 
though its derivative is not zero. As the value of 0o decreases the maximum and the minimum 
approach each other. 

(iii) For a /2b > (f>a > 0, the two extrema (that collapse into a saddle point at ~ a /2b) 
start separating again as 00 is reduced. The maximum, located now at 0m — 0o returns to 
at 00 = 0, while the minimum, 0„i = a/b — 0o, increases as 0o is decreased. 

Observe that in the interval [0, a/b] the effective potential has the same form at a/2b + A 
and at a/26 — A. Note also that in that interval the potential has two minima, one with zero 
derivative, and another at the origin with, in general, non-vanishing derivative, which is due 
to the presence of an absorbing barrier at = 0. It is also possible to decide the fraction 
of total density located at each minima; supposing in first approximation that the potential 
is locally minimized at all the points simultaneously, and calling p the fraction at = and 
I ~ p the complementary probability in such a way that 0o = Op + (1 — p)4'mi we obtain 
p — {a — 260o)/(a — 60o). In particular the mixture is symmetric at 0o — a/{Zb). 

We have performed some numerical simulations to further study the proposed model. For 
that we consider an initial condition with average field density 0o and small fluctuations around 
that value. In figure 1 we show a long time typical density profile in d = 1 {L = 100) for two 
different initial densities in the unstable regime: 0o = 0.20 and 0o — 0.26. The two phases 
are clearly asymmetric; the phase at = is free of fluctuations. At 0o — 0.26 ( i.e. slightly 
below the critical point of instability, a/26 = 0.3 for the chosen parameter values) the fraction 
of sites in the nonvanishing region is larger than in the other case, but in any case, the walls 
are quite soft. 

We have verified the presence of only one characteristic length is the scaling regime. For 
that we have computed the structure function at different times, and rescaled it using R{t)] a 
good collapse of the curves is observed (see Fig. 2). 

For the time evolution of the averaged cluster size we get, in d = 1: R{t) ~ ^("±0.02)^ with 
a — 0.30,0.28 and 0.26, for initial concentrations 0o = 0.05,0.1 and 0.15 respectively (Fig. 3). 

Deciding whether the assymptotic behavior is governed by an exponent 1/3 or 1/4 is a 
difficult task from both theory and numerics. From the theoretical point of view it has been 
argued that when the bulk dynamics is suppressed an exponent 1/4 can be obtained |13] , while 
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Fig. 1. - Density profile in the one-dimensional case for (f)Q = 0.20 (lower curve) and 0o = 0.26 (upper 
curve) respectively at t ^ 10^. The parameter values are taken to be: a = 0.6, = 1, and D = 0.01. 
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Fig. 2. - Rescaled structure function in d = 1 for times from t = 4000 to t — 24000. 



time dependent exponents, interpolating between 1/3 and 1/4 (reflecting crossover effects) are 
obtained for concentration dependent diffusion coefficient [y_3| . Observe that in our model the 
dynamics ceases in the interior of absorbing clusters but not in the other (active) phase. More 
systematic simulations and further theoretical analysis aiming to sort out this point will be 
presented elsewhere. 

In d = 2 a similar qualitative behavior is obtained. In Fig. 4 we show the evolution of 
the field variable as a function of time for an initial density 0o = 0.05. Observe that first 
absorbing (black) regions are created in a sea of active sites (white) (at time t = the whole 
lattice would be white) . Afterwards these islands percolate through the system and the active 
regions become isolated. 

All the measured effects in both d = 1 and d = 2, are essentially deterministic, and 
correspond therefore to the T — fixed point of the conserved RFT. Therefore, the main 
ingredient characterizing this phase separation in systems with an absorbing state, is the 
physical constraint at (j) — 0, and not the noise itself. 

Our system may exhibit a phase transition from the ordered behavior to a disordered phase 
as the noise amplitude is increased. Due to the presence of a noise of a different nature we 
expect this transition to be in a universality class other than that of model B. This issue will 
be analyzed in a forthcoming paper. 

Summing up, we have shown that deterministically monostable systems, and, in particular. 
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Fig. 3. - Log-log plot of the averaged mean cluster radius as a function of time for d = 1, and different 
initial densities. 
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Fig. 4. - Time evolution of the density in a 100 x 100 lattice for an initial average concentration 
(j)o = 0.05 (with small fluctuations around that mean value). Black points represent sites in the 
absorbing state, while active regions (i.e. above a certain tiny threshold) are represented in black. 
First the initial condition in unstable and absorbing regions are created it = 50). Then these regions 
coalesce it = 80) and percolate through the system, isolating regions of active sites (£ = 100). The 
size of active regions shrinks down for late times (t — 150). 

systems with absorbing barriers can exhibit the phenomenon of phase separation, in a form 
that shows many common features with the standard picture of phase separation, but has 
also some essential differences. These come mainly from the fact that the two minima of the 
effective potential are intrinsically of different nature in our model. 

Among other many possible applications of the general phenomena we have just described 
let us cite the following: phase separations of charged particles in a conductor (observe that 
these situation has an absorbing phase: zones with no charge are fluctuation free), or diffusion 
of particles in an absorbing substrate. We expect the general scenario studied here to describe 
some realistic model similar to the ones cited previously. 

Finally, coming back to our toy model in economics, the wealth remains homogeneously 
distributed among individuals when the total wealth in the society is larger than a critical 
value, (determined by the government tax policy), but when it is smaller than the critical 
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threshold the society becomes separated in two groups: rich and poor, which separate from 
each other. 

It is a pleasure to acknowledge useful discussions with J. M. Sancho J. Marro, P. Garrido, 
L. Pietronero, and R. Toral. This work has been partially supported by the European Union 
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